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PREFACE 


The  motivation  of  ysis,  performed  in  support  of  the  Air 

Force’s  efforts  in  thet.  n’ssile  defense  (TMD),  was  to  bring  in¬ 
creased  attention  to  the  development  and  fielding  of  capabilities 
to  hold  at  risk  enemy  theater  ballistic  missiles  (TBMs)  and  their 
transporter-erector-launchers  (TPLs)  through  anacks  on  the  TELs 
either  before  or  after  TBM  launch. 

The  report  should  be  of  interest  to  a.nlysts  in  ;he  Tf.fD  community 
and  to  others  interested  generally  in  the  application  of  probabilistic 
and  Monte  Carlo  methods  to  military  problems. 

The  study  was  sponsored  by  the  Air  Combat  Command  Arrpament 
Requirements  Division  (ACC/DRA)  and  took  place  within  the  Force 
Modernization  and  Employment  Program  of  Project  AIR  FORCE,  the 
Air  Force’s  federally  funded  research  and  development  center  at 
RAND. 

PROJECT  AIR  FORCE 

Project  AIR  FORCE,  a  division  of  RAND,  is  the  Air  Force  federally 
funded  research  and  development  center  (FFRDC)  for  studies  and 
analyses.  It  provides  the  Air  Force  with  independent  analyses  of 
policy  alternatives  affecting  the  development,  employment,  combat 
readiness,  and  support  of  current  and  future  aerospace  forces.  Re¬ 
search  is  being  performed  in  three  programs;  Strategy,  Doctrine,  and 
Force  Structure;  Force  Modernization  and  Employment:  and  Re¬ 
source  Management  and  System  Acquisition. 
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Project  AIR  FORCE  is  operated  under  Contract  F49620-91-C-0003 
between  the  Air  Force  and  RAND. 

Brent  Bradley  is  Vice  President  and  Director  of  Project  AIR  FORCE. 
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SUMMARY 


Under  the  assumption  that  future  opponents  will  choose  to  acquire 
inventories  of  theater  ballistic  missiles  (TBMs)  substantially  larger 
than  the  number  of  their  transporter-erector- launchers  (TELs),'  this 
report  supports  the  view  that  counterforce  operations  against  the 
TELs  can  play  a  role  in  sharply  reducing  the  overall  size  of  a  prospec¬ 
tive  threat.  For  initial  inventory  ratios  of  10  TBMs  per  TEL  (a  ratio 
not  inconsistent  with  past  practice  by  many  countries),  reductions  of 
approximately  80  percent  in  missiles  launched  are  possible  with 
probabilities  of  successful  postlaunch  TEL  kill  of  about  0.5.  Even  for 
probabilities  of  TEL  kill  of  only  0.2,  reductions  of  50  percent  are  pos¬ 
sible.  Combined  prelaunch  and  postlaunch  counterforce  attacks  act 
synergistically,  enhancing  the  overall  effectiveness.  The  report  does 
not  discuss  the  circumstances  in  which  these  attacks  can  achieve  any 
specific  degree  of  effectiveness.  History  suggests  that  claims  for 
significant  counterforce  capabilities  should  be  viewed  with  skepti¬ 
cism.  While  we  concur  with  this  observation,  there  may  exist  a 
significant  motivation  for  striving  to  make  this  capability  at  least 
modestly  effective.  We  believe  that  this  is  clearly  possible  for 
postlaunch  counterbattery  operations. 

It  is  useful  to  note  that  postlaunch  counterforce  attacks  (often  called 
counterbattery  attacks)  do  little  initially  to  affect  TBM  launches.  As 


*This  assumption  is  reasonable  for  current  threats  because  the  cost  of  the  TFl.s  is 
significantly  greater  than  the  cost  of  the  TBMs.  In  the  future,  as  missiles  get  smarter 
and  TELs  therefore  can  become  less  complex,  the  TELs  may  become  single-shot 
“throw-aways."  For  the  foreseeable  future,  however,  we  would  expect  to  see  many 
reloads  per  TEL,  although  not  necessarily  many  tons,  as  was  the  case  in  Iraq. 
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such,  counterbattery  capabilities  are  properly  considered  to  be  ad¬ 
juncts  to  various  active  defense  measures  for  defeating  TBM  threats. 
They  aid  these  active  defensive  measures  mainly  by  lessening  the 
magnitude  of  the  attack  (both  in  terms  of  the  total  number  of  TBMs 
that  can  be  launched  and  their  launch  rate),  thus  potentially  lower¬ 
ing  the  total  investment  in  active  defenses. 

The  report  also  provides  the  reader  with  a  set  of  equations  (and  their 
derivations)  to  calculate  additional  outcomes.  For  those  interested 
in  such  calculations,  an  appendix  lists  computer  programs  suitable 
for  this  purpose. 
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_ Chapter  One 

INTRODUCTION 


The  purpose  of  this  report  is  to  bring  increased  attention  to  the  de¬ 
velopment  and  fielding  of  capabilities  to  hold  at  risk  enemy  theater 
ballistic  missiles  (TBMs)  and  their  transporter-erector-launchers 
(TELs)  through  attacks  on  the  TELs  either  before  or  after  TBM 
launch.  The  "conventional  wisdom"  seems  to  be  that  postlaunch  (or 
"counterbattery”)  attacks  on  the  TELs  are  feasible  (because  of  the 
exploitable  missile  launch  signature),  but  not  very  useful.  After  all, 
the  process  starts  with  the  successful  launch  of  the  threat  missile,  and 
empty  TELs  have  little  intrinsic  value  in  themselves.  On  the  other 
hand,  successful  prelaunch  attacks  against  the  TBMs/TELs  would  be 
highly  useful  (killing  the  TBMs  before  they  were  launched,  thereby 
avoiding  a  host  of  end-game  concerns),  but  are  regarded  as  techni¬ 
cally  and  operationally  far  more  challenging  and  thus  have  an  ex¬ 
tremely  low  probability  of  success.  Thus,  conventional  wisdom  fur¬ 
ther  suggests  that  efforts  in  both  areas  are  likely  to  have  low  payoff. 
This  study  challenges  that  assessment.  We  do  not  dispute  that  these 
concepts  are  challenging,  although  we  do  agree  that  counterbattery 
attacks  are  inherently  easier  to  achieve  than  prelaunch  kills.  But  we 
do  believe  that  even  modest  capabilities  could  be  very  useful.  If 
those  capabilities  could  be  fielded  with  relatively  little  expense  (as  we 
believe  is  the  case  by  adapting  existing  systems  such  as  the  F-15E 
Strike  Eagle),  then  our  analysis  supports  the  view  that  they  could  be 
"best  buys.”’ 


’This  statement  assumes  that  the  United  States  has  control  of  the  airspace  over 
potential  TBM  operating  areas.  Assuming  this  to  be  the  case,  the  additional  option  of 
boost-phase  kill  also  exists.  The  synergy  between  prelaunch,  boost-phase,  and 
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In  recent  work  by  David  Vaughan  and  Richard  Mesic,  the  value  of  a 
counterbattery  attack  was  represented  by  the  number  of  residual 
TBMs  that  could  no  longer  be  launched  if  TELs  were  located  and  de¬ 
stroyed  immediately  after  their  associated  TBMs  were  launched.^ 
Figure  1  shows  the  results  of  a  calculation  that  was  used  to  describe 
this  result.  In  this  simplistic  calculation  it  was  assumed  that  there 
were  an  infinite  number  of  potential  reloads,  so  that  the  average 
number  of  successful  launches  could  be  derived  fi'om  the  geometric 
probability  distribution  (i.e.,  the  expected  number  is  1 /postlaunch 
TEL  kill  probability).  In  Figure  1  we  simply  truncated  these  curves  at 


RANDMn«».| 


Postlaunch  TEL  kill  probability 

Figure  1 — Counterbattery  Attacks  Can  Ground  Reload  Missiles 


postlaunch  attacks  against  the  missile  and  the  TEL  helps  make  all  three  more  cost- 
effective. 

^Unpublished  work  by  David  Vaughan  and  Richard  Mesic  on  investment  strategy  for 
theater  missile  defense;  unpublished  briefing  by  Tim  Nafi'  (U.S.  Army  Space  and 
Stratepc  Defense  Command  ISSDC)),  “TMD  Attack  Operations  and  Active  Defense 
Functional/Quantitative  Relationships." 
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the  assumed  total  TBM  inventory  of  1000.  The  Naff  briefing  made 
the  point  that  even  relatively  small  single-engagement  kill  probabili¬ 
ties  could  substantially  reduce  the  total  size  of  the  attack  by  killing  all 
the  TELs  before  all  the  TBMs  could  be  launched.  While  a  full  com¬ 
plement  of  active  terminal-area  defenses  would  still  be  required  to 
handle  the  TBMs  that  were  launched,  the  total  attack  against  those 
defenses  would  be  reduced.* 

The  Naff  briefing  supported  the  feasibility  of  obtaining  a  competent 
counterbattery  capability.  Not  only  would  the  exhaust  plume  of  the 
in-flight  TBM  provide  an  excellent  cue  for  alerting  the  searching  sen¬ 
sors  to  the  presence  of  a  TEL,  but  backtracking  the  missile’s  trajec¬ 
tory  to  locate  a  launch  point  on  the  ground  would  provide  a  relatively 
small  footprint  area  within  which  the  TEL  could  be  located.  Various 
practical  operational  concepts  exist  for  prosecuting  the  cue  and  the 
backtrack  into  a  successful  attack  against  the  empty  TEL,  including 
the  employment  of  F-15Es  with  suitable  software  modifications  to 
the  APG-70  radar.^ 

Thus,  Figure  1  supports  arguments  in  favor  of  obtaining  a  robust 
counterbattery  capability.  However,  it  begs  some  additional  ques¬ 
tions: 

•  Figure  1  is  based  on  an  approximation  (i.e.,  it  assumes  an  infinite 
supply  of  TBMs).  Would  more  exact  (hence,  more  complex) 
analyses  such  as  we  develop  here  for  more  realistic  inventory 
sizes  and  operational  constraints  (e.g.,  TEL  and  TBM  clustering 
and  finite  numbers  of  TBMs)  show  a  reduced  utility  for  counter¬ 
battery  attacks? 

•  How  would  the  results  change  if  prelaunch  counterforce  attacks 
were  included?  Although  the  actual  capability  of  prelaunch  at- 


*The  reduction  in  threat  size  is  sensitive  to  the  ratio  of  the  attacker’s  TBMs  to  TELs. 
The  greater  the  ratio,  the  higher  is  the  utility  of  postlaunch  counterbattery  attacks. 
While  the  ratio  of  TBMs  to  TELs  is  obviously  writhin  the  control  of  the  opponent, 
historically  countries  have  chosen  many  more  TBMs  than  TELs.  Ratios  of  10  or  larger 
are  the  rule. 

^See,  for  example,  Mesic  and  Vaughan,  “A  Concept  for  Near-Tem  Action  on  TMD 
Counterforce  Capabilities  Development  and  Demonstration,”  May  1993,  and  Mesic, 
Vaughan,  and  Shaver,  "Theater  Missile  Defense  (TMD)  Pillar  Balancing;  A  Strawman 
Approach  and  Initial  Observations,”  April  1993  (project  memoranda  to  the  Air  Force). 
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tacks  might  be  signincantiy  less  than  postlaunch,  TELs  would  be 
killed  in  both. 

•  Given  the  random  nature  of  the  outcomes,  with  what  confidence 
could  we  assert  a  particular  outcome,  given  that  we  knew  (or 
were  willing  to  guess  at)  the  likelihood  of  successful  location, 
identification,  and  kill?  If  counterbattery  attacks  rest  on  killing 
all  the  TELs,  the  tails  in  the  distribution  of  outcomes  M  be 
extremely  important. 

We  address  these  questions  in  an  analytic  context.  The  foJi,)vving 
chapters  will  first  cover  counterbattery  attacks  when  only  a  single 
TEL  is  involved,  then  extend  the  analysis  to  cover  multiple  TELs  de¬ 
ployed  as  part  of  a  cluster  of  TBMs  and  TELs,  and  then  further  extend 
the  analysis  to  include  prelaunch  counterforce  attacks.  A  few  obser¬ 
vations  follow.  Appendix  A  contains  supporting  mathematics  that 
should  help  anybody  wishing  to  pursue  these  subjects.  Appendix  B 
lists  various  computer  programs  (written  in  QuickBASIC)  that  were 
developed  to  support  the  figures  shown  in  this  report. 
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Chapter  Two 

THE  SIMPLE  MATHEMATICS  FOR  COUNTERBATTERY 
ATTACK  OUTCOMES  ASSUMING  A  SINGLE  TEL 


Assume  that  a  single  TEL  has  access  to  N  TBMs — one  on  the  TEL  it¬ 
self  and  an  additional  N  - 1  back  at  the  resupply  site.  Further  assume 
that  if  the  TEL  is  killed,  whatever  TBMs  remain  at  the  resupply  site 
can  no  longer  be  launched  (later  we  will  assume  that  other  TELs  can 
access  these  TBMs,  an  assumption  that  will  affect  our  outcomes). 
Then  the  probability  that  a  TEL  successfully  launches  exactly  K  TBMs 
is  simply  the  probability  that  it  survives  the  counterbattery  attack 
after  the  first  K  - 1  launches  but  is  killed  after  the  Kth  launch,  that  is, 


1<K<N  (1) 


and 


P(N)  =  (l-Pk)'''‘  K  =  N  (2) 

where  is  the  probability  that  the  TEL  is  killed  after  any  particular 
launch. 

Using  Eqs.  (1)  and  (2),  we  can  derive  the  following  formula  for  the 
expected  number  of  TBMs  launched,  E(K): 


E,(K) 


(3) 
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where  the  numerical  subscript  after  E  identifies  the  size  of  the  TEL 
inventory  (in  this  case,  one).  The  derivation  of  Eq.  (3)  can  be  found 
in  Appendix  A.  For  very  large  N  or  very  small  TEL,  survival  probabil¬ 
ity  (for  small  values  of  (1  -  Pt)^),  E,(K)  becomes  simply 

the  approximated  formula  used  in  the  earlier  Vaughan-Mesic  work. 
In  addition,  if  we  let  P|^  approach  zero,  then  Ei(K)  approaches  N; 
that  is,  if  the  TELs  are  never  killed,  they  launch  the  full  inventory  of 
TBMs.  Figure  2  plots  the  expected  size  of  the  TBM  inventory 
launched,  Ei(K),  as  a  function  of  for  various  TBM  inventory  sizes. 
Shown  for  comparison  is  the  expected  number  of  TBMs  launched, 
assuming  an  infinite  supply  of  TBMs,  so  that  the  nature  of  the  errors 
in  that  simple  approximation  can  be  seen.  As  the  counterbattery  kill 
potential  falls,  the  divergence  between  the  finite  and  infinite  TBM 


.01  .1  1 
Countertattery  probability  of  kill 


Figure  2— Ability  of  Counterbattery  Attacks  to  Limit  the  Expected  Number 
of  TBMs  Launched,  Assuming  a  Single  TEL 


Counterbattery  Attack  Outcomes  Assuming  a  Single  1  tL  7 
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stockpiles  grows,  with  the  finite  stockpile  outcome  asymptoting  to  its 
inventory  size  as  approaches  zero. 

Two  points  are  worth  noting: 

•  For  TBM  inventories  per  TEL  greater  than  about  5,  counterbat¬ 
tery  Pj,s  as  modest  as  0.2  can  dramatically  reduce  the  total  frac¬ 
tion  of  the  TBM  inventory  that  can  be  launched.  Other  work 
suggests  that  P^s  of  0.2  or  higher  (perhaps  much  higher)  may  be 
achievable  at  only  modest  cost  and  risk.' 

•  High  counterbattery  kill  potential  can  limit  the  number  of  TBMs 
launched  even  if  the  inventory  per  TEL  is  relatively  small.  Of 
course,  if  there  is  only  one  TBM  per  TEL,  then  counterbattery 
attacks  have  no  analytic  effect  (although  the  effect  on  the  morale 
and  motivation  of  the  doomed  launch  crews  might  be  signifi¬ 
cant). 

Because  of  the  random  nature  of  these  outcomes,  it  is  also  of  interest 
to  look  at  the  variance  of  these  outcomes.  Figure  3  plots  the  statisti¬ 
cal  variance  (the  square  of  the  standard  deviation)  of  these  outcomes 
as  a  function  of  P|j  for  several  values  of  N.  Note  that  the  variance  can 
be  relatively  large  at  intermediate  values  for  P|,  (peaking  where  Pk  = 

I /the  number  of  TELs),  leading  to  concerns  that  simply  using 
expected  values  for  counterbattery  attack  outcomes  in  multilayer 
attack  calculations  may  not  properly  capture  the  consequences  of 
unlucky  occurrences.  For  example,  for  P^  =0.1  and  10  TBMs  per 
TEL,  the  expected  number  of  launches  (from  Figure  2)  is  about  6.5, 
but  the  standard  deviation  (the  square  root  of  the  curves  in  Figure  3) 
is  about  4 — a  significant  number  compared  with  the  expected  value. 


'See.  for  example,  Richard  Mesic.  “Extended  Counterforce  Option.s  for  Coping  uith 
Tactical  Ballistic  Missiles,"  in  Paul  Davis  (ed.).  New  Challenges  for  Defense  Planning, 
RAND,  Santa  Monica,  CA,  1994. 
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Counterbattery  probability  of  kill 


Figure  3— Variance  of  Counterbattery  Attack  Outcomes  for 
2-10  TBMs  and  a  Single  TEL 


Chapter  Three 

EXTENDING  THE  OUTCOMES  TO  CASES  WHERE 
MULTIPLE  TELs  HAVE  ACCESS  TO  A  COMMON 

CACHE OFTBMs 


The  simple  model  of  a  single  TEL  with  TBM  reloads  that  it  alone  can 
access  and  launch  is  probably  not  a  realistic  portrayal  of  how  most 
enemies  would  operate.  More  likely  is  a  model  of  a  cache  of  TBMs 
and  a  number  of  TELs  that  have  access  to  that  cache.  Unless  all  TELs 
with  access  to  that  cache  are  killed,  the  remaining  TBMs  in  the  cadie 
can  still  be  launched. 

The  mathematics  associated  with  calculating  expected  outcomes  in 
cases  where  TBMs  and  TELs  are  clustered  into  caches  of  multiple 
TELs  (T)  and  TBMs  (N)  is  provided  in  Appendix  A.  No  simple  closed- 
form  solutions  equivalent  to  Eq.  (3)  exist,  and  solving  the  equations 
must  be  done  numerically  on  computers.  Three  approaches  were 
used.  The  first  approach  yields  a  set  of  algebraic  equations  that, 
while  not  reducible  to  a  simple  closed  form,  can  be  easily  calculated 
by  numerical  means.  This  approach  provides  exact  answers  for  both 
the  expected  values  of  the  outcomes  and  the  associated  distribution 
functions  for  the  number  of  TBMs  actually  launched  for  any  specific 
set  of  inputs.  It  also  provides  the  reader  with  an  understanding  of 
how  the  primary  input  variables  affect  the  outcome  without  per¬ 
forming  extensive  calculations.  However,  the  numerical  calculations 
can  be  time-consuming  if  the  number  of  TELs  in  a  cluster  is  large 
(greater  than  5),  and  thus  is  only  useful  for  small  T. 

A  second  approach  is  to  computationally  solve  a  matrix  calculation 
patterned  after  Markov.  This  approach  also  produces  exact  answers 
for  both  the  expected  values  of  the  outcomes  and  the  associated  dis¬ 
tribution  functions  for  the  number  of  TBMs  actually  launched  for 
any  specific  set  of  inputs.  It  is  not  computationally  limited  to  small 
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numbers  of  TELs;  it  can  rapidly  calculate  outcomes  for  essentially 
any  threat  size. 

The  third  and  last  approach  uses  a  Monte  Carlo  computer  simulation 
to  yield  estimates  (rather  than  exact  calculations)  of  expected  out¬ 
comes  and  data  on  outcome  distributions:  it  also  works  for  large 
threat  sizes.  If  the  number  of  trials  is  sufficiently  large,  the  uncer¬ 
tainties  in  outcome  estimates  resulting  from  the  stochastic  character 
of  the  approach  can  be  minimized. 

Table  1  summarizes  these  characteristics.  Appendix  B  contains  the 
computer  codes  developed  for  each  of  these  approaches. 

The  next  several  figures  address  the  question  of  the  importance  of 
assuming  that  TELs  and  TBMs  are  clustered  into  sets  wherein  every 
TEL  has  access  to  and  can  launch  every  TBM.  It  is  possible  for  there 
to  be  multiple  clusters,  in  which  case  we  assume  that  a  TEL  from  one 
cluster  cannot  launch  TBMs  that  belong  to  another  (if  they  could, 
then  we  would  simply  increase  the  size  of  the  cluster).  An  extreme 
case  is  if  all  TELs  have  access  to  all  TBMs.  For  practical  reasons,  this 
seems  unlikely.  However,  the  other  extreme — each  TEL  can  access 
only  its  specifically  assigned  reload  TBMs  (i.e.,  all  clusters  consist  of 
single  TELs,  as  mathematically  portrayed  in  Chapter  Two) — seems 
equally  unlikely. 

Figure  4  starts  with  the  assumption  that  the  total  number  of  TBMs  is 
ten  times  the  total  number  of  TELs.  It  shows  how  clustering  TELs 
diminishes  the  effectiveness  of  postlaunch  counterbattery  attacki, 
normalizing  the  outcome  to  the  TBM/TEL  ratio. 

Table  1 

Mathematical  Approaches  to  Calculating  Counterbattery 
Attack  Outcomes 


Calculation 

Method 

Calculated  Outcomes 

Expected  Values  Distributions 

Computer 

Limitations 

Algebraic 

Exact 

Exact 

Requires  small 
number  (T)  of 
TELs  for  cluster 

Markov  matrix 

Exact 

Exact 

None 

Monte  Carlo 

Approximate 

Approximate 

None 

Extending  the  Outcomes  to  Multiple  TELs  with  Common  TBM  Cache  1 1 
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Counterbattery  probability  of  kill 


Figure  4 — Expected  Number  of  TBMs  Launched  (Normalized  to  Total 
Number  of  TELs)  as  a  Function  of  the  Number  of  TELs  in  a  Cluster 


Several  observations  can  be  drawn  from  Figure  4: 

•  If  the  counterbattery  kill  probability  is  above  about  0.2,  then  the 
clustering  of  TELs  does  not  significantly  affect  the  outcome. 

•  Clustering  ratios  above  about  5  have  diminishing  returns. 

•  Regardless  of  clustering,  counterbattery  kill  probabilities  above 
about  0.2  can  still  sharply  diminish  the  total  number  of  TBMs 
that  can  be  launched. 

Another  important  variable  is  the  ratio  of  total  TBMs  to  TELs. 

Assuming  a  cluster  of  100  TBMs,  Figure  5  shows  the  expected 

number  of  these  TBMs  that  could  be  launched  as  a  function  of  the 

number  of  TELs  added  to  the  cluster.  We  see  the  following: 

•  An  obvious  and  effective  counter  to  counterbattery  attacks  is  to 
buy  additional  TELs.  The  sensitivity  is  substantial,  but  even 
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Counterbattery  probability  of  kill 


Figure  5 — Sensitivity  of  Counterbattery  Attack  Effectiveness 
to  TEL  Proliferation 


buying  a  fairly  large  number  of  TELs  does  not  totally  obviate 
counterbattery  effectiveness. 

•  If  P|( 's  of  0.5  or  higher  are  achievable,  even  a  large  number  of 
TELs  cannot  guarantee  that  the  majority  of  TBMs  in  the  inven¬ 
tory  can  be  successfully  launched. 

It  is  impossible,  of  course,  to  judge  from  these  calculations  whether 
TEL  proliferation  constitutes  a  cost-effective  response  to  a  counter¬ 
battery  threat.  Past  experience  with  regard  to  TBM  deployments 
suggests  that  countries  have  decided  the  issue  by  buying  a  large 
number  of  TBMs  per  TEL.  No  doubt  this  decision  was  influenced  in 
part  by  the  perception  that  future  opponents  could  not  mount 
successful  counterbattery  attacks.  But  it  almost  certainly  was  also 
affected  by  peacetime  operating  costs  and  other  constraints. 
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Chapter  Four 

AN  EXAMPLE  ANALYSIS,  ASSUMING  5  TEU  AND  50 

TBMs  IN  A  CLUSTER 


The  following  example  demonstrates  how  an  analysis  of  counter¬ 
battery  operations  might  be  carried  out 

Assume  that  the  threat  consists  of  an  enemy  who  has  500  TBMs  and 
50  TELs,  divided  into  10  clusters.  Each  cluster  is  assumed  to  operate 
autonomously.  The  questions  we  address  are  (1)  on  the  average, 
how  much  would  counterbatterv  attacks  reduce  the  overall  si^e  of 
the  threat?  (2)  with  what  confidence  can  we  state  these  outcomes? 
and  (3)  how  sensitive  are  the  answers  to  assumptions  about  how  the 
enemy  operates  and  the  probability  that  counterbattery  attacks  will 
succeed? 

Figure  6  shows  •;he  expected  number  of  TBMs  that  would  be 
launched  from  a  specific  cluster  as  a  function  of  counierbattery  ef¬ 
fectiveness:  it  simply  repeats  one  of  the  curves  in  both  Figures  2  and 
4.  It  could  have  been  calculated  by  any  of  the  three  approaches,  al¬ 
though  the  first  two  give  an  exact  answer.  For  comparison,  we  have 
included  the  infinite  TBM  case.  As  noted,  counterbattery  kill  prob¬ 
abilities  above  about  0.2  result  in  a  strong  likelihood  that  the  TELs 
would  be  killed  prior  to  TBM  exhaustion,  thus  sharply  curtailing  the 
total  number  of  TBMs  that  would  be  launched. 

Also  shown  in  Figure  6  are  the  10  percent  and  90  percent  outcomes. 
The  lower  boundary  is  defined  by  specifying  that  90  percent  of  the 
outcomes  exceed  it,  whereas  only  10  percent  exceed  the  upper 
boundary.  The  curves  show  that  outcomes  can  vary  widely. 
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Counterbattery  probability  of  kill 


Figure  6 — Expected  Values  Versus  Upper  and  Lower  Boundaries, 
SOTBMsandSTELs 


If  we  had  constructed  similar  curves  for  the  full  inventory  of  500 
TBMs,  the  boundaries  would  have  become  narrower.  We  have  not 
calculated  the  results,  but  they  can  be  approximated  by  adjusting  the 
differences  by 

VlO, 

where  10  comes  from  the  ratio  of  the  total  inventory  to  the  number  of 
TBMs  in  the  cluster.' 

Figure  7  is  shows  the  distribution  of  outcomes  for  this  case  another 
way.  This  curve  has  been  calculated  for  =  0.2.  Note  that  the 


'Assuming  ample.TELs,  the  variance  of  the  results  about  the  mean,  normalized  by  the 
mean,  go  as  1  /  V n ,  where  n  is  the  total  number  of  TBMs.  As  TEL  survival  comes  into 
play,  the  variance  will  change.  However,  assuming  a  factor  of  10  in  increased  numbers 
of  both  TBMs  and  TELs,  the  normalized  variance  of  the  larger  inventory  should  be 
about  1/ VIO. 
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RAHDUfi469-7 


Figure  7 — Fraction  of  Cases  Where  Specific  Number  of  TBMs 
Is  Launched,  Pk  =  0.2 


maximum  likelihood  outcome  is  about  20  TBMs,  although  the  mean 
is  just  under  25.  The  "jump”  or  spike  at  the  end  corresponds  to  the 
fraction  of  trials  where  all  of  the  TBMs  are  launched.  This  spike  is 
large  whenever  the  counterbattery  is  sufficiently  small  that  mul¬ 
tiple  TELs  survive  to  TBM  exhaustion.  Note  that  the  spread  in  the 
outcome  is  substantial,  as  we  could  have  anticipated  from  Figure  6. 

For  comparison,  we  also  ran  the  Monte  Carlo  model  for  this  set  of 
inputs.  Table  2  shows  the  outcomes  of  the  run.  The  numbers  in  the 
histograms  are  the  occurrences  where  the  given  number  of  TBM 
launches  or  surviving  TELs  was  the  result  (the  number  of  TBMs 
launched  can  run  from  a  minimum  of  5  to  a  maximum  of  50;  the 
surviving  TELs  from  0  to  5).  For  example,  out  of  the  500  trials  run  for 
this  case,  in  482  all  the  TELs  were  killed,  and  in  20  trials  exactly  26 
TBMs  were  launched  (read  ‘6’  across  the  top  and  ‘20’  down  the  side 
to  find  the  number  of  cases  where  26  TBMs  were  launched).  For  this 
specific  run,  the  calculated  average  number  of  TBMs  fired  is  slightly 
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Table  2 

Output  from  Monte  Carlo  Model  (Counterbattery  Attacks)  I 


Input  Data 

Total  TELs  in  cluster  =  5 
Total  TBMs  in  cluster  =  50 
Pkl-  0 
Pk2  =  0.2 

Number  of  trials  =  500 
Output  Data 

Average  number  of  TBMs  fired  =  26.136 

Standard  deviation  =  10.109 

Standard  deviation  of  mean  estimate  =  0.452 


Histogram  of  Number  of  TBMs  Fired 
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Histogram  of  Number  of  TELs  Surviving 
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higher  than  the  theoretical  expected  outcome  (in  fact,  given  the 
standard  deviation  of  the  mean  of  the  estimate,  this  case  is  a  sub¬ 
stantial  outlier).  Figure  8  plots  the  TBM  data  in  Table  2  on  top  of  the 
data  in  Figure  7. 

What  can  we  derive  from  this?  In  some  plausible  scenarios,  a  reduc¬ 
tion  in  both  the  total  size  of  the  threat  and  in  the  salvo  sizes  possible 
at  any  particular  time  could  gready  enhance  TMD  effectiveness;  in 
other  scenarios  the  effects  might  be  small.  Certainly,  counterbattery 
attacks  of  this  character  do  little  to  deny  a  threat  based  on  weapons 
of  mass  destruction,  where  denial  of  any  successful  TBM  penetra¬ 
tions  is  the  dominant  criterion.  To  deal  with  such  threats,  we  need  to 
consider  both  postlaunch  active  defense  options  (e.g..  Theater  High- 
Altitude  Area  Defense  (THAAD)  or  boost-phase  defenses)  and 
prelaunch  counterforce  options. 
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Figure  8 — Monte  Carlo  Outcomes  Superimposed  on  Figure  7,  Pjj  =  0.2 


_ Chapter  Five 

EXTENDING  THE  CALCULATIONS  TO  INCLUDE 
PRELAUNCH  COUNTERFORCE  ATTACKS 


The  above  calculations  presumed  that  the  only  threat  to  the  TELs 
and  TBMs  came  from  attacks  against  empty  TELs  after  the  TBM  had 
been  launched.  It  is  obvious  that  the  same  aircraft  that  might  exe¬ 
cute  the  postlaunch  counterbattery  attack  will  be  in  a  position  to 
carry  out  prelaunch  counterforce  attacks  if  the  TBM /TEL  can  be 
found  prior  to  launch.  The  mathematics  for  prelaunch  counterforce 
attacks  can  be  found  in  Appendix  A.  We  will  not  repeat  the  math 
here.  However,  for  comparison  with  Eq.  (3),  we  offer  the  general 
formula  for  the  expected  number  of  TBMs  launched  assuming  a  sin¬ 
gle  TEL. 


Ei(K) 


(5) 


where  we  have  designated  the  prelaunch  single-engagement  kill 
probability  as  ‘  to  differentiate  it  from  the  engagement  kill  prob¬ 
ability  associated  with  postlaunch  attacks  (we  will  designate  post¬ 
launch  single-engagement  kill  probabilities  by  A  comparison 
with  Eq.  (3)  shows  that  the  two  formulas  are  identical  with  the  excep¬ 
tion  of  the  factor 


in  the  numerator.  Note  that  when  the  probability  of  TEL  survival 
over  N  launches  is  very  small  (i.e.,  when 
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is  small),  Eq.  (5)  shows  that,  as  is  intuitively  obvious,  one  less  launch 
would  be  expected  than  in  the  postlaunch  case  in  Eq.  (3). 

Figure  9  is  equivalent  to  Figure  2,  except  it  assumes  that  there  are 
only  prelaunch  counterforce  attacks  against  a  single  TEL.  Note  that 
these  curves  are  slightly  more  favorable  to  the  attacker,  because  a  kill 
results  in  denying  a  TBM  launch  as  well  as  killing  the  TEL. 

Figure  9  depicts  the  no  TEL  clustering  case.  Figure  10,  the  equivalent 
of  Figure  4,  shows  the  impact  of  TBM /TEL  clustering  on  the 
prelaunch  counterforce  outcome.  Again,  the  similarity  to  Figure  4  is 
not  surprising,  with  counterforce  attacks  of  equal  kill  potential 
yielding  slightly  more  favorable  results  for  the  attacker.  However,  we 
note  that  prelaunch  counterforce  attacks  are  not  blessed  with  the 
same  certainty  of  a  cueing  signal  as  are  the  postlaunch  counter- 
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Prelaunch  countarforce  probability  of  kill 


Figure  9 — ^Ability  of  Prelaunch  Counterforce  Attacks  to  Limit  the  Expected 
Number  of  TBMs  Launched,  Assuming  a  Singie  TEL 
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Prelaunch  counterfbrce  probability  of  kill 


Figure  10— Effect  of  Clustering  on  Prelaunch  Counterforce 
Attack  Outcomes 


battery  attacks.  In  the  latter  case,  the  bright  infrared  (IR)  signal 
associated  with  TBM  launch  is  virtually  impossible  to  mask  or 
otherwise  conceal.  Finding  TELs  without  such  a  cue  proved  to  be 
almost  impossible  during  the  Gulf  War,  and  it  is  not  certain  that 
significant  progress  has  been  made  since  then  to  reverse  this 
unfavorable  outcome. 

Combined  prelaunch  and  postlaunch  counterforce  attacks  would 
obviously  improve  the  outcomes.  The  attacker  would  get  two  oppor¬ 
tunities  to  kill  the  TEL  during  a  launch  sequence,  effectively  increas¬ 
ing  the  kill  probability  as  well  as  killing  TBMs  if  the  prelaunch  attack 
succeeds.  The  expected  number  of  TBMs  launched  for  a  single 
TEL — the  equivalent  equation  to  Eqs.  (3)  and  (5) — is  simply 
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where  the  prescripts  1  and  2  denote  prelaunch  and  postlaunch. 

Figure  11  displays  outcomes  from  combined  pre-  and  postlaunch  at¬ 
tacks,  where  it  is  assumed  that  the  cluster  consists  of  100  TBMs  and 
10  TELs.  We  have  parameterized  the  postlaunch  attack  probability, 
plotting  curves  for  different  prelaunch  capabilities. 

The  results  in  Figure  1 1  show  a  striking  synergy  between  prelaunch 
and  postlaunch  attack  effectiveness.  Even  modest  prelaunch  capa¬ 
bilities  (e.g.,  PkS  in  the  neighborhood  of  0.2)  can  signiticantly  aug¬ 
ment  the  effectiveness  of  postlaunch  counterbattery  attacks,  reduc¬ 
ing  the  number  of  TBMs  that  can  be  launched  by  50  percent  or  more. 

Clearly,  the  countermeasure  of  adding  more  TELs  would  lessen  the 
effectiveness  of  the  combined  attacks.  In  general,  however,  com¬ 
bined  attacks  can  sharply  curtail  the  total  launch  potential  of  a  TBM- 
armed  enemy.  In  addition,  counterforce  attacks,  particularly 
prelaunch  attacks,  will  aid  theater  missile  defenses  by  lessening  the 
number  of  TBMs  launched  during  any  launch  period  as  well  as 
overall.  Counterforce  attacks  in  combination  with  boost-phase 


NAIIOUn46»-l> 


Postlaunch  counterbattery  probability  of  kill 


Figure  1 1 — Combined  Effectiveness  of  Prelaunch  and  Postlaunch 
Counterforce  Attacks:  100  TBMs  and  10  TELs  in  Cluster 
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attacks,  attacks  that  could  be  carried  out  by  the  same  aircraft,  would 
make  any  conceivable  enemy  countermeasures  to  most  theater 
missile  defense  architectures  much  less  effective.  Thus,  counterforce 
attacks  can  directly  (as  well  as  indirectly)  support  overall  theater 
missile  defense  architectures,  improving  the  effectiveness  of  other 
components  of  the  defense. 

We  return  to  the  earlier  example  of  50  TBMs  and  5  TELs.  Table  3  re¬ 
peats  Table  2,  with  an  overlay  of  a  prelaunch  counterforce  attack 
effectiveness  of  0.1. 
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Tables 

Output  from  Monte  Carlo  Model  (Counteibattery  Atuudts)  II 


Input  Data 

Total  TELs  in  cluster  =  5 
Total  TBMs  in  cluster  =  50 
Pj.l=0.10 
P^2  =  0.20 

Number  of  trials  =  500 
Output  Data 

Average  number  of  TBMs  fired  =  16.378 

Standard  deviation  =  7.185 

Standard  deviation  of  mean  estimate  =  0.321 


Histogram  of  Number  of  TBMs  Fired 
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Chapter  Six 

OBSERVATIONS 


Representative  calculations  of  the  potential  implications  of  post¬ 
launch  and  prelaunch  counterforce  attack  operations  against  TELs 
and  their  associated  TBMs  indicate  that  counterforce  operations  can 
substantially  reduce  the  overall  threat  of  TBM  attacks.  It  can  achieve 
this  effect  even  if  the  probabilities  of  success  for  individual  attacks 
are  modest.  For  example,  with  initial  inventory  ratios  of  10  TBMs  per 
TEL  (a  ratio  not  inconsistent  with  past  practice  by  many  countries), 
reductions  of  approximately  80  percent  in  missiles  launched  are  pos¬ 
sible  with  probabilities  of  successful  postlaunch  TEL  kill  of  about  0.5. 
Even  for  probabilities  of  TEL  kill  of  only  0.2,  reductions  of  50  percent 
are  possible.  Combined  prelaunch  and  postlaunch  counterforce  at¬ 
tacks  act  synergistically,  enhancing  the  overall  effectiveness.  We 
have  not,  however,  discussed  the  circumstances  in  which  these  at¬ 
tacks  can  achieve  any  degree  of  effectiveness.  History  suggests  that 
claims  for  significant  counterforce  capabilities  should  be  viewed  with 
skepticism.  While  we  concur  with  this  observation,  this  discussion 
suggests  that  there  exists  a  significant  motivation  for  striving  to  make 
this  capability  at  least  modestly  effective.  We  believe  that  this  is  pos¬ 
sible  for  postlaunch  counterbattery  operations. 
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Appendix  A 

THE  MATHEMATICS  OF  COUNTERFORCE  ATTACKS 


Appendix  A  describes  the  mathematics  behind  estimating  the  poten¬ 
tial  benefits  of  counterforce  attacks  (both  prelaunch  attacks  against 
enemy  theater  ballistic  missiles  (TBMs)  and  their  transporter- 
erector-launchers  (TELs),  and  postlaunch  attacks  against  just  the 
TELs).*  The  discussion  will  take  place  in  three  sections.  The  first  will 
develop  the  formulas  for  calculating  the  expected  number  of  TBMs 
that  could  be  launched  as  a  function  of  (1)  the  likely  success  of 
individual  counterforce  attacks  (Pk),  (2)  the  number  of  TELs  in  a 
TEL/TBM  cluster,  and  (3)  the  number  of  available  TBMs  in  the  same 
cluster.  These  formulas  usually  cannot  be  stated  in  closed  form  and 
require  solution  by  computer.  Because  of  their  nature,  this  limits  the 
magnitude  of  the  numbers  of  TELs  and  TBMs  that  can  be  calculated 
with  these  formulas. 

The  second  section  will  describe  a  Markov  matrix  approach  to  solv¬ 
ing  the  same  problem.  This  approach  also  requires  computer  solu¬ 
tion.  but  alleviates  some  of  the  constraints  on  the  numbers  of  TELs 
and  TBMs  that  can  be  included  in  the  calculation. 

The  third  section  will  describe  a  Monte  Carlo  computer  model  that 
also  solves  this  problem.  This  approach  has  the  benefit  of  relatively 
rapid  calculation  of  outcomes  for  large  numbers  of  TBMs  and  TELs, 
with  the  constraint  that  the  outcomes  are  inherently  random. 
Outputs  from  a  number  of  runs  will  be  described. 


^Common  usage  is  to  call  prelaunch  attacks  counterforce,  and  postlaunch  attacks 
counterbattery.  We  will  generally  follow  this  usage. 
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SECTION  ONE:  SOLVING  THE  PROBLEM  BY  ALGEBRAIC 
METHODS 

Part  One:  The  Mathematics  of  Postlaunch  Attacks  Against 
the  TELs  Only 

A.  Assumptions.  We  start  with  the  following  assumptions. 

1.  The  TBMs  and  TELs  are  organizationally  formed  into  clusters.  A 
cluster  is  a  unit  of  TBMs  and  TELs  that  operates  autonomously. 
For  purposes  of  reloading,  every  TEL  in  the  cluster  is  assumed  to 
have  access  to  every  TBM  in  the  cluster.  Once  all  TBMs  in  the 
cluster  are  used,  the  cluster  is  empty  (i.e.,  it  is  not  resupplied). 
Also,  if  all  TELs  belonging  to  the  cluster  are  killed,  any  remaining 
missiles  in  the  cluster  are  lost. 

2.  The  number  of  TELs  in  the  cluster  is  T. 

3.  The  number  of  TBMs  in  the  cluster  is  N.  N  is  larger  than  T 
(reloads  are  assumed  to  exist). 

4.  Each  TEL  is  loaded  with  a  TBM  when  it  leaves  its  hide/resupply 
site.  The  TELs  are  not  subjected  to  attrition  on  their  way  to  the 
launch  site  and  are  assumed  to  successfully  launch  the  TBM  at  the 
site.  Counterbattery  attacks  are  assumed  to  occur  while  the  TEL  is 
returning  to  its  temporary  hide  or  reload /resupply  site. 

5.  All  TELs  are  treated  the  same.  Each  survives  the  attempted 
counterbattery  attack  (after  every  launch)  with  probability  (1  -  P|(). 

These  assumptions  reflect  the  operational  concept  for  counter¬ 
battery  attacks — that  the  approximate  location  of  the  TEL  is  de¬ 
termined  by  backtracking  the  TBM  while  it  is  in  powered  flight  (i.e., 
while  it  has  a  very  bright,  readily  detectable  signal  that  would  be 
difflcult  to  conceal  or  mask),  and  an  airborne  platform  suitably 
armed  flies  to  that  location,  tries  to  detect  and  identify  the  (now 
empty)  TEL,  and  prosecutes  the  attack.  The  mathematics  are 
generally  applicable  to  alternative  concepts  of  operations,  but  we  will 
use  this  one  for  ease  of  description. 

B.  The  Mathematics  for  a  Single  TEL  per  Cluster.  Let  P(K)  =  proba¬ 
bility  that  exactly  K  TBMs  are  successfully  launched.  This  can  hap¬ 
pen  only  if  the  TEL  successfully  launches  its  TBM,  recovers  and 
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reloads  (K  -  1)  times  and  is  subsequently  killed  after  the  Kth  launch 
The  only  exception  to  this  Mould  be  the  case  where  the  Kth  launch  is 
equal  to  N.  Thus  we  nddres"  two  separate  cases,  one  where  K  <  N 
and  one  where  K  =  N. 

ForK<N: 

P(K)  =  (l-Pkf'‘*  Pi,  (A.1) 

where 

is  the  probability  that  the  TEL  survives  the  first  K  -  1  recoveries  from 
the  launch  site  and  Pi,  is  the  probability  that  it  gets  killed  on  the  last. 
Note  that  K  runs  from  1  to  N,  reflecting  the  fact  that  the  first  launch  is 
not  contested  (all  counterbattery  attacks  occur  after  the  launch)  and 
the  TEL  cannot  launch  more  missiles  than  exist  in  the  cluster. 

ForK  =  N: 

P(K)  =  (l-Pk/^'‘’  (A.2) 

In  other  words,  to  launch  all  N  TBMs,  the  TEL  needs  only  to  survive 
the  first  (N  -  1)  counterbattery  attacks.  Whether  it  survives  the 
counterbattery  attack  after  the  last  launch  makes  no  difference  to  the 
number  of  TBMs  launched  (i.e.,  the  cluster  is  empty). 

Expected  value  for  1  TEL  and  N  TBMs:  Let  E,{K)  =  expected  number 
of  TBMs  that  are  successfully  launched.  Ei(K)  is  given  by  the  follow¬ 
ing  sum: 

E,{K)=  XK  P{K)=  i'k  Pk  (l-Pkf"' 

K=I  K=1 

N-1 

{'\.3) 


+  N  (l-Pk) 
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where  Eqs  (A.l)  and  (A.2)  have  been  substituted  into  the  first  sum¬ 
mation.  To  solve  this  equation  in  closed  form,  recall  the  following 
series  expansions: 

- —  =  1-e2x-e  3x^  +  4x^  +  ...  =  yix‘“'  (A.4) 

(1-xf  ifr 

^  - 1 1 X I  x^  I  I  =  y  x'“' 


Also  note  the  following  identity: 


Xf(i)=if(i)-Xf(i)  (A.6) 

i=I  i=l  i=N 

Using  Eqs  (A.4),  (A.5),  and  (A.6), 

j,i  x  (l-x)'  '  =  x£i  (l-x)''‘  -xX>  (l-x)’’*  (A.7) 

i  =  l  i=l  i=N 

To  solve  these  two  new  summations,  note  that  the  first  summation  is 
simply  the  expansion  shown  in  Eq.  (A.4).  The  second  summation  has 
a  similar  character,  but  the  series  does  not  start  at  1.  To  solve  it, 
change  variables  in  the  last  summation,  letting  j  =  i  -f  1  -  N.  This 
yields 

i=i 

(A.8) 

=  i-x  (l-xf''£(N-l-Ej)(l-x)''' 

^  j=i 


(l-(l-x)) 
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Expanding  the  remaining  sum  into  two  new  sums,  we  get 

i  (A.9) 


where  we  have  used  Eqs.  {A.5)  and  (A.4)  on  the  two  sums,  respec¬ 
tively.  If  we  now  substitute  the  pertinent  variables,  that  is,  =  x 
and  N  =  n,  into  equation  (A.9),  and  substitute  this  into  Eq.  (A.3),  we 
obtain 


E,(K) .  N(i- p,)”-' -iXi- ?,)“■' 


(A.  10) 


which  simplifies  to  our  final  form: 


E,{K)  =  — 


{A.11) 


Note  that  if  N  is  large,  the  second  term  in  the  numerator  becomes 
nearly  zero,  yielding  the  result: 


LimE,(K)  =  —  (A.12) 

This  simple  expression  for  the  expectation  can  easily  be  derived  from 
probability  theory.  For  our  purposes,  however,  Eq.  (A.11)  is  of  more 
interest,  because  we  generally  believe  that  TBM  stockpiles  will  be  fi¬ 
nite  and  not  out  of  proportion  to  the  number  of  available  TELs, 

In  addition  to  the  expected  value,  it  is  useful  to  calculate  the  variance 
(or  standard  deviation).  If  Var(K)  is  the  variance  of  these  outcomes, 
then 


I 
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Var{K)=  |;K2p(K)-(E(K)f 

K  =  1 


(A.  13) 


After  substituting  Eqs.  {A.1),  {A.2),  and  (A.10)  into  Eq.  (A.13),  we  find 
that  no  simple  closed -form  solution  exists.  Thus, 


Var(K)  =  N2(i-  P|cf 

K  =  l 


1-(1-P,f 


(A.  14) 


Note  that 


Lim  Var(K)  =  0 
f)c->o 


and 


Lim  Var(K)  =  0, 


which  is  as  it  should  be. 

C.  The  Mathematics  for  2  TELs  per  Cluster.  The  mathematics  in¬ 
volving  multiple  TELs  per  cluster  adds  the  complication  that  there 
are  multiple  ways  to  obtain  exactly  K  successful  launches.  As  above, 
let  P(K)  equal  the  probability  that  exactly  K  TBM  launches  are  made. 
Then 

For2  <K<N; 


K-l 

P(K)  =  |Prob(TEL  #  Haunches  exactly  i  TBMs)  | 

i=l 


{Prob{TEL  #  2  launches  K  -  i  TBMs) } 


(A.15) 


i 

1 


I 


I 
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which,  from  above,  we  know  leads  to 

p(K) = X  {(i  -  Pk  r'  Pk }  •  {(i  -  pk  Pk } 

„  ,  (A.  16) 

=i{i-Pkrv=(K-i)Pk^(i-Pkr 

i=I 

Note  that  K  starts  at  2,  reflecting  the  assumption  that  both  TELs 
launch  their  first  TBMs  without  being  contested. 

ForK  =  N; 


P(N)  =  1- X'PCK) 


K=2 


(A.  17) 


which,  after  substituting  Eq.  (A.16)  into  {A.17),  becomes 


P(N)  =  1  -  X  -  Pk  f  (A.18) 

K=2 

which,  after  a  suitable  change  of  variables  (j  =  K  - 1),  simplifies  to 


p(N)=i-p/xi(^-Pkr’ 

1=1 


(A.  19) 


Before  proceeding  to  further  refine  Eq.  (A.19),  let  us  back  up  to  Eqs. 
(A.16)  and  (A.18)  and  write  the  expression  for  the  expected  number 
of  TBMs  launched.  As  before. 


E2(K)  =  X  K  P(K)  =  N  -  NPk^  X(*^  -  Pkf 

K=2  K=2 

+  Pk'lK(K-l)(l-Pkf'' 

K=2 


(A.20) 
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which  reduces  to 


E2(K)  =  N  -  Pk^  2;'(N  -  k)(K  -  l)(l  -  Pkf (A.21) 

K=2 

This  form  for  the  expected  value  will  reappear  later  when  the  sum¬ 
mation  terms  cannot  be  readily  simplified. 

Equation  (A.  19)  is  similar  to  Eq.  (A.7),  and  can  be  simplified  in  the 
same  manner.  Expanding  into  two  summations  and  changing  vari¬ 
ables  in  the  second  (to  get  the  summation  from  one  to  infinity),  we 
obtain 


N-2 

i=t 


1-1 


=ij(i-Pkr-i(j-N-2Xi-Pkr"- 

j=i  j=i 


(A.22) 


The  first  term  in  Eq.  (A.22)  is  simply 


1 


To  solve  the  second  summation,  separate  it  into  two  new 
summations  as  follows: 


io.N-2)(i-Pkr^^*-’=(i-p,r{ij(i-Pkr' 

i=i  ij=i 

+i(N-2)(i-p.r'! 


(A.23) 


(A.24) 


Substituting  Eq.  (A.24)  back  into  Eq.  (A.22)  and  that  sum  back  into 
Eq.  (A.19),  we  obtain 
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P(N)  =  (l -  Picf '"{l  +  Pk(N  -  2)}  (A.25) 

Thus,  by  suitable  algebraic  manipulation  we  have  managed  to  obtain 
closed-form  solutions  for  the  probabilities  of  exactly  K  TBMs  being 
launched  for  all  values  of  K. 

Expected  value  for  2  TELs  and  N  TBMs:  Substituting  Eq.  (A.25)  into 
the  formula  for  the  expected  value  of  the  number  of  TBMs  launched 
yields 

E^IK)  =  X  K  •  P{K)  =  n(1-  '"{l+  P^N  -  2)} 

IC=2  ^ 

(A.26) 

+  XK(K-l)P,^(l-P,f'" 

K=2 

Alas,  the  authors  are  not  aware  of  a  closed-form  solution  for  this 
summation.  However  the  equation  is  easily  calculated  by  numerical 
means,  and  the  authors  have  built  a  simple  BASIC  program  for  this 
purpose  (see  Appendix  B).  Also  note  that  Eq.  (A.26)  is  not  much 
simpler  to  calculate  than  Eq.  (A.21). 

D.  The  Mathematics  for  3  TELs  per  Cluster. 

For  3  S  K  <  N;  The  above  approach  of  deriving  the  probability  of 
exactly  K  launches  will  continue  to  be  followed,  but  with  3  TELs  it  is 
necessary  to  introduce  another  variable.  Thus,  let  I,,  I2,  and  I3  be  the 
number  of  TBMs  launched  by  TELs  1,  2,  and  3,  respectively. 
Obviously,  for  any  K, 

1,  +  I2  +  I3  =  K  (A.27) 

Thus, 
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P(K)  =  I Prob(TEL#llaunchs exactly  I,  TBMs) 
ii=i 


K-i-h 


[  Prob(TEL  #  2  launches  exactly  I2  TBMs)  ■ 
Prob(TEL  #  3  launches  I3  TBMs)]J 


l2  =  l 


And  recognizing  as  before  that 

Prob(TEL  #1  launches  exactly  li  TBMs)  =  Pit(l  -  Pr)  * 

Prob(TEL  #2  launches  exactly  I2  TBMs)  =  P|j(l  -  P^ )  ^ 

Prob(TEL  #3  launches  exactly  I3  TBMs)  =  P^  (l  -  P^  )'^  ^ 

then 


P(K)  =  I  Pk{1-  Pk)'‘""lk{l-  Pky^'‘Pk(l-  Pk)^'""^-’ 

Ii=i  I2=1  (A.28) 

Pulling  out  the  appropriate  terms  from  the  summations,  we  obtain 

„  ,K-2  K-h-!  , 

p(K)=p,'(i-Pk) ■  i  i=Pk(i-Pk)  i(K-i-ii) 


,K-2  K-li-! 

I  i 

II  =  l  l2  =  l 

Expanding  the  summation  yields 


(A.29) 


^In  calculating  the  sum  or  a  series  of  consecutive  integers,  it  is  usef'il  to  recall  the 
“trick”  attributed  to  Gauss  when  he  was  in  the  first  grade.  Add  the  sum  to  itself  with 


I 


I 


I 

I 

1 
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For  K  =  N: 

By  derinition,  and  using  Eq.  (A.30)  for  P(K)  where  K  <  N. 


N-l 

P(N)=  1-  XP(K) 

K=3 

(A.31) 

K=3  * 

Expected  Value  for  3  TELs  and  N  TBMs:  Substituting  into  the 
expression  for  the  expected  value  of  number  of  TBMs  launched, 

EjCK)  =  X  K  P(K) 

K=3 


N-t 


K=3 

N-l 


3/.  „,k-3(K-l](K-2) 


K=3 


(A.32) 


3N  (K-1)(K-2), 

=  N-Pk  X(N-K)i - - ^1-Pk 

K=3  2 


,K-3 


The  authors  do  not  know  of  a  way  to  solve  this  series  in  a  simpler 
form.  However,  computer  codes  have  been  built  to  evaluate  it  (see 
Apendix  B). 


the  series  reversed:  that  is,  the  first  number  is  added  to  the  last,  the  second  to  the  next 
to  last,  etc.  The  resulting  sums  are  all  constant,  and  there  are  exactly  n  of  them,  where 
n  is  the  number  of  integers  in  the  series.  Therefore,  the  new  sum  is  just  the  product  of 
the  number  of  integers  times  the  sum  of  the  first  and  last  term  in  the  series.  To  get  the 
answer,  simply  divide  by  two,  reflecting  the  fact  that  you  have  added  the  series  twice. 
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E.  The  Mathematics  of  4  TELs  per  Complex.  The  approach  taken  is 
identical  to  that  in  D.  above — solve  first  for  the  probability  of  exactly 
K  launch  successes  (K  less  than  or  equal  to  the  total  number  of  TBMs 
in  the  complex),  and  then  plug  the  solution  into  the  expected-value 
equation.  Thus, 

For  4  <  K  <  N : 


K-3, 


P(K)=  5^{Prob(TEL 

ii=i 


#1  launches  exactly  Ij  TBMs) 


K-l-Ii 

^  Prob(TEL  #  2  launches  exactly  I2  TBMs) 

l2  =  l 


K-UI1-I2. 

^  I  Prob(TEL  #  3  launches  exactly  I3  TBMs) 

13=1 


Prob(TEL#  4  launches  exactly  I4  TBMs)|| 


K-3  ,  ,K-2-Ii  ,  , 

ll  =  l  l2  =  I 

13=1 


.  ^  .K-3  K-2-I1  K-l-Ii-h 

=>v'(i-p.rz  i  t 

ll  =  l  l2  =  l  13=1 

.  „  .K-3  K-2-I1 

Il  =  l  l2  =  l 
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) 


which  yields 


P{K)=  lv"(l-Pk 


(K-l-Ii)(K-2-U) 

2 


(A.33) 


ForK  =  N: 


N-l 

P(N)  =  1-  XW) 

K=4 

K=4  !l  =  l  ^ 

Expected  value  for  4  TELs  and  N  TBMs: 

As  above, 


E,(K)  =  N-X(N-K)Pk'(l-I^,)' 

K=4 


K-4 


k^3(k-1-I,XK-2-I,) 

2 

ii=i  ^ 


(A.35) 


There  is  no  way  to  signihcantly  simplify  this  expression.  It  is,  how¬ 
ever,  readily  amenable  to  numerical  calculation. 

F.  Generalization  of  Mathematics  to  Larger  Numbers  of  TELs  per 
Cluster.  The  generalization  should  be  obvious.  The  reader  can 
confirm  the  following  equations: 
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ForK<N: 


P(K)=P/(1-Pkf''' 


[k-(T-1)|k-(T-2)-I,Ik-(T-3)-1,-12) 

I  s  s 

II  >2  >3 


lT-1 


(A.36) 


For  K  =  N; 


p(N)  =  i-  |;p/(i-Pkf^ 


-(T-1)Ik-(T-2)-I,] 

I  I 


K=T  I)  l2 

[K-(t-3)-I,-I2]  [K-1-1,-  -It-2] 

S  I 


>3 


'T-1 


(A.37) 


If,  for  each  notation,  we  define  the  summations  as  Sb(K,T),  where  the 
subscript  b  denotes  that  this  sum  applies  for  counterbattery  calcula¬ 
tions,  then 


[k-(T-1)][K-(T-2)-I||k-(T-3)-Ii-12)  [K-1-I,-  -IT-2J 

Sb(K.T)=  X  X  X  X  (A.38) 

l|  I2  '3  <T-1 


and 


Et(K)  =  N  -  X(N  -  K)PJ(1-  Pkf  ■■'’Sb(K,T)  (A.39) 

K=T 


Equation  (A.39)  is  the  general  solution  to  this  expected-value  prob¬ 
lem.  The  equation  is  amenable  to  computer  calculation,  although 
the  time  required  for  its  calculation  grows  substantially  as  T  and  N 
increase.  For  practical  purposes,  calculations  for  T  >  10  are  probably 
not  best  done  by  this  equation. 


( 
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Part  Two:  Extending  the  Mathematics  for  Prelaunch  Attacks 
Against  Both  the  TBMs  and  TELs 

A.  Assumptions.  The  first  three  assumptions  in  Part  One  are  the 
same.  However,  because  the  attack  can  kill  both  the  TELs  and  the 
TBM,  changes  to  assumptions  4  and  5  are  needed. 

4.  Each  TEL  possesses  a  TBM.  The  TELs  are  subjected  to  attrition  on 
their  way  to  the  launch  site.  A  successful  prelaunch  counterforce 
attack  will  kill  both  the  TBM  and  the  TEL. 

5.  All  TBM/TELs  are  treated  the  same.  Each  survives  the  attempted 
counterforce  attack  (before  each  launch)  with  probability  (1  -  Pi,). 

B,  The  Mathematics  for  a  Single  TEL  per  Cluster.  Much  of  the  for¬ 
malism  is  identical  to  that  used  above.  The  primary  changes  concern 
the  range  for  the  summations  (they  start  with  zero  instead  of  one.  as 
it  is  possible  that  no  TBMs  will  be  launched).  We  again  start  by 
looking  at  the  probability  that  exactly  K  TBMs  are  launched. 

For  0  <  K  <  N: 


P(K)  =  {1-Ptf  Pk 


(A.40) 


where  the  TEL  successfully  survives  the  first  K  sorties  out  of  its  hide 
site,  launching  K  TBMs,  but  is  discovered  and  killed  on  the  K  +  1  sor¬ 
tie. 

ForK=N: 


P(N)  =  {l-P,f 


(A.41) 


where  the  TEL  survives  all  N  sorties,  launching  all  N  TBMs. 
Expected-value  equation: 

E,(K)=  XKP(K)=N(l-P|,f +  x'K(l-Pk)'Pk  (A.42) 

K=0  K=0 
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The  summation  at  the  end  of  Eq.  (A.42)  can  be  solved  in  closed  form 
by  the  techniques  used  aboxe^).  Solving  this  equation  yields  the  fol¬ 
lowing  result: 


E,(K)  =  |i-l 
Pk  ) 


(A.43) 

Comparing  this  with  Eq.  (A.ll)  shows  the  addition  of  a  -1  term  to 


Pk 


reflecting  the  fact  that  killing  the  TEL  results  in  loss  of  one  TBM  as 
well. 

C.  The  Mathematics  for  2  TELs.  Much  of  the  above  applies  for  2 
TELs  as  well.  However,  there  is  at  least  one  important  difference. 
Because  a  TBM  is  killed  every  time  that  a  TEL  is  killed,  it  is  important 
to  consider  separately  the  case  where  exactly  1  (out  of  2)  TELs  sur¬ 
vives  to  TBM  exhaustion. 

As  before,  for  0  <  K  <  (N  - 1) 

P(K)  =  t(l-P,)‘p,{l-P,f'‘p, 

i=0 


=  {K.1)P,^(1-P,f 

For  K  N  -  1,  one  of  the  TELs  must  be  killed,  and  one  must  have 
survived  to  TBM  exhaustion.  Thus, 

p(N-i)= x(i-p,yp,(i-Pkf" 


•^Extract  a  (l  - 1\ )  from  the  summation. 
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where  the  kill  of  the  first  TEL  can  occur  anytime,  including  up  k.  the 
time  when  the  TEL  is  carrying  the  last  TBM  to  the  launch  stand. 
Thus, 


P(N-1)=  NPk(l  -  P,,)'" 


tA.45) 


Finally,  tor  all  TBMs  to  be  launched,  all  TELs  must  survive.  Thus, 


(A.46) 


Using  Eqs.  (A. 44),  (A.45),  and  (A.46),  the  expected  value  can  be  writ¬ 
ten  as 


E.(K)  =  N(1  -  -c  N(N  -  l)l\(l  -  ' 

.,  (A. 47) 

+  +  l)l\'(l  -  l\f 

K=0 

This  equation  cannot  be  significantly  simplified. 

D.  The  Mathematics  for  3  TELs.  The  equations  for  2  TEl.s  point  the 
way,  but  now  we  must  consider  as  separate  cases  the  situation  where 
exactly  1  and  2  TELs  survive  to  TBM  exhaustion. 

If  all  TELs  are  killed  before  TBM  exhaustion,  then 

For  K  <  N  -  2: 


pw = 1(1-  p^yp/iV-  po’p.(i  -  p.r^'p. 

i=0  j=0 


For  K  =  N  -  2  (2  TELs  killed.  1  survives  to  TBM  exhaustion): 
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i 


+  I)(K  H 

K-l 


(A.52) 


E.  The  Mathematics  for  N  TELs.  The  symmetry  in  the  above  equa¬ 
tions  should  now  be  evident.  Using  the  notation  from  the  first  part  of 
this  section,  we  can  now  write  a  general  expression  for  the  various 
probabilities  and  the  expected  value. 
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ForO  <  K  <  N-T: 


P(K)  =  (1-P,)’'PJX  X 

I|=0  l2=0 

-*T-3  .  . 

I  {K  +  l-I,-l2 - 1^.2) 

'T-2 

Define  a  new  variable  for  the  summation:  i.e., 

K  K-h 

Sf(K.T)=£  X  •  • 

tl=0  l2=0 

-lT-3  , 

I  {K  +  l-I,-l2 - 1^.2) 

It -2 


K  K-h 


(A.53) 


(A.54) 


then 


P{K)  =  p/(l-P,,fSf(K,T) 


(A.55) 


For  K  =  N  -  (T  -  1),  where  all  but  one  TEL  is  killed  before  TBM  ex¬ 
haustion, 


P(K)  =  pJ''(l-PkfSf(K,T) 


(A.56) 


ForK  =  N-  (T-  2), 

P(N-(T-2))=  pJ-"(l-P,f''''-^’Sf{K,T-l)  (a.57) 

And,  in  general,  for  N  -  T  <  K  <  N, 

P(N  -  (T -,•))=  pJ-‘(l-P,f"'^‘''Sf(K,T-j  +  l)  (A.58) 


where  I  <  j  <  T  -  1 
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ForK  =  N. 


P(N)  =  (l-P,f 


Part  Three:  The  Mathematics  of  Combined  Prelaunch  and 
Postlaunch  Counterforce  Attacks 

Because  there  are  simpler  ways  than  the  algebraic  approach  to  solve 
the  combined  prelaunch  counterforce  and  postlaunch  counterbat¬ 
tery  attacks,  we  will  consider  only  the  single  TEL  case. 

Let  P,  be  the  probability  of  successful  prelaunch  counterforce  attack 
against  the  TEL,  and  P2  the  probability  of  successful  counterbattery 
attack  against  the  TEL  after  it  has  launched  the  TBM.  Let 
Qi=(l-P[)  and  Q2=(l-P2)-  In  addition,  define  Q3=QiQ2. 
Then,  for  K  <  N, 

P(K).{l-P,)'‘(l-P,)'‘-'jp,+(l-P,)P,} 

where  the  only  way  that  exactly  K  TBMs  can  be  launched  is  for  the 
TEL  to  survive  K  prelaunch  counterforce  engagements  and  K  -  1 
postlaunch  engagements,  only  to  be  killed  (the  terms  in  the  brackets) 
after  the  kth  TBM  launch  by  a  counterbattery  attack  or  by  a 
prelaunch  engagement  on  the  TEL’s  way  to  launch  its  K  -t-  1  TBM. 
This  equation  can  be  rewritten  into  the  following  form: 

P(K)  =  Q,{Q,Q2f"‘[l-Q,Q2] 


or 


P(K)  =  Q,Q3''-'(1-Q3)  (A.61) 

For  K  =  N,  the  formula  is  simply 


•  P{N)  = 


(A.62) 


I 


The  Mathematics  of  Counterforce  Attacks  47 


As  above,  we  can  derive  the  expected  number  of  TBMs  launched  by 
suitable  algebra,  resulting  in  the  following  equation: 


E,(K) 


Qi(i-Q3'') 

(1-Q3) 


(A.63) 


A  quick  comparison  between  this  equation  and  those  for  counterbat¬ 
tery  and  counterforce  shows  that  this  formula  is  identical  to  them  if 
either  kill  probability  is  set  to  zero. 


Section  Two:  The  Markov  Matrix  Approach  for  Obtaining 
Outcomes  for  Counterforce  and  Counterbattery  Attacks^ 

A  different,  and  in  many  ways  more  satisfying,  approach  to  solving 
the  problem  is  through  the  approach  which  might  roughly  be 
equated  to  a  Markov  matrix  model.  To  describe  this  approach,  we 
make  the  following  assumptions. 

1.  At  the  start  T  TELs  and  N  TBMs  are  assumed  to  exist  and  to  be 
colocated  at  a  common  hide  site  (for  analytic  ease,  we  will  assume 
that  the  so-called  operational  hide  site  and  the  resupply  site  are 
the  same).  We  assume  that  N  >  T,  that  is  there  is  more  than  one 
TBM  per  TEL. 

2.  Each  step  in  the  process  starts  by  assuming  that  a  single  TEL, 
outfitted  with  a  TBM,  leaves  its  hide  site,  committed  to  launching 
the  TBM.  On  route  to  the  launch  site,  the  TEL  and  its  TBM  are 
subjected  to  a  counterforce  attack.  If  the  attack  is  successful,  both 
the  TEL  and  the  TBM  are  killed.  If  the  attack  fails,  we  assume  that 
the  TEL  successfully  reaches  the  launch  site  and  launches  the 
TBM.  After  the  TBM  is  launched,  the  TEL  attempts  to  return  to  its 
hide  site.  We  assume  that  the  TEL  is  again  subject  to  attack,  and 
successfully  returns  to  its  hide  only  if  it  survives  this  attack  as  well. 
If  the  TEL  reaches  its  hide  site  successfully,  it  is  available  to  be 
reloaded  with  a  new  TBM  and  reused. 


■'The  author  owes  a  debt  of  gratitude  to  David  Vaughan  for  suggesting  this  approach. 
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3.  The  TEL  is  not  attacked  except  on  its  way  to  and  from  the  launch 
site. 

4.  If  an  individual  TEL  is  killed,  a  replacement  TEL  (if  the  initial 
stockpile  of  T  TELs  has  not  been  exhausted)  repeats  the  above 
process. 

5.  The  process  terminates  when  all  T  TELs  have  been  killed  or  all  the 
TBMs  launched. 

Part  One:  The  Postlaunch  Counterbattery  Attack  Case  Only 

For  ease  of  presentation,  we  start  with  the  postlaunch-counterbat- 
tery-attack-only  case.  Let  t(i,x)  be  the  probability  of  exactly  i  TELs 
surviving  after  x  attempts  to  launch  a  TBM.  Assume  for  simplicity 
thatT  =  3.  and  let  P,  =  (l  -  P,j).  Then 

't(3,x)1  TPs  0  0  Ol  rt{3,x-l)' 

t{2,x)  ^  Pk  Ps  0  0  t(2,x-l) 

t(l,x)  0  Pk  Ps  0  ^  t(l,x-l)  (A.64) 

t{0,x)J  [o  0  Pk  ij  [t(0,x-l) 

with  the  following  initial  conditions: 

t(3,0) 
t(2,0) 
t(1.0) 
t(0,0) 

Define  E(x)  as  the  expected  number  of  TBMs  launched,  given  x  at¬ 
tempts.  Thus, 

T 

E(x)  =  ^  t(i,  X  -  1)  +  E(x  - 1) 

i=l 


(A.66) 
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Alternatively,  E(x)  can  be  written  as 

E(x)  =  (l-t(0,x-l))  +  E(x-l)  (A.67) 

More  generally,  assume  that  A(i,j)  is  the  probability  matrix  that  trans¬ 
forms  t(j,x-  1)  to  t(i,x).  Then 

t(i,x)=A(i,j)  t(j,x-l)  (A.68) 

Because  the  matrix  A  is  relatively  sparse,  it  is  computationally  sim¬ 
pler  to  pose  the  problem  in  the  form  of  the  following  set  of  difference 
equations: 


For  i  =  T 

t(T,x)=  P,t(T,x-l) 

(A.69) 

For  0  <  i  <  T 

t(i,  x)  =  Pjtd,  X  - 1)  -1-  P|(t(i  -1- 1,  X) 

(A.70) 

For  i  =  0 

t(0,x)=  Pkta,x-l)-^t(0,x-l) 

(A.71) 

Finally,  x  runs  from  I  to  N. 

Part  Two:  Extending  the  Markov  Matrix  Model  to  Include 
Prelaunch  Counterforce  Attacks 

The  above  approach  is  readily  extended  to  include  prelaunch  coun¬ 
terforce  attacks.  Assuming  that  prelaunch  and  postlaunch  attacks 
are  successful  with  probabilities  *  P|j  and  ^  P|( ,  respectively,  then 


For  i  =  T 


t(T,x) 


t(T.x-l) 


(A.72) 


For  0  <  i  <  T 


t(i,x)  =  (l-'Pk)(l-2p,,)t(i,x-l)9-|'P|,  -t^l-'P|^)2Pk  t(i  +  l,x-l) 


(A.73) 
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Fori  =  0  t(0.x)  =  t(0,x-l)-E|^P|, -E^l-lPkj^Pkjta.x-l)  (A.74) 

The  coefficients  of  t{i,x)  are  the  elements  of  the  Markov  matrix.  The 
bracketed  term  can  be  rewritten  into  a  somewhat  simpler  form;  i.e., 

{'P.*(i-'P.)'P>)-(i-Q.Q,)-(i-Q,) 

where  Qj  and  Q2  are  the  TEL  survival  probabilities  against 
prelaunch  and  postlaunch  counterforce  attacks,  respectively,  and 
Qj  =  Q1Q2.  Substituting  Q,,  Q2,  and  Qj  into  Eqs.  (A.72),  (A.73),  and 
{A.74),  we  obtain  the  following  compact  form: 

Fori  =  T  t(T,x)  =  Q3t(T,x-l)  (A.75) 

For  0  <  i  <  T  t(i,  x)  =  Qstfi,  x  -  1)  -i-  (l  -  Q3  )t(i  + 1,  x  -  1)  (A.76) 


Fori  =  0  t(0,x)  =  t(0,x-l)-E{l-Q3)t(l,x-l)  (A.77) 

The  number  of  TBMs  launched,  E(x),  takes  a  slightly  different  form: 

E(x)  =  Q,  •  (1  - 1(0, x-l)]+  E(x -  1)  (A.78) 


where  a  TBM  is  successfully  launched  only  if  the  TEL  survives  the 
pre-launch  counterforce  attack.  Also  note  that  the  probability  that 
exactly  x  TBMs  are  launched  is  just  the  probability  that  after  x  launch 
attempts  the  last  TEL  is  killed.  If  we  set  f(x)  as  this  probability,  then 

f(x)  =  t(0,  X)  - 1(0,  X  -  1)  (A.79) 


Part  Three:  Solving  the  Markov  Matrix 
It  is  interesting  to  note  that 

t(i,x)=  AMD)  t(j,0)  (A-80) 
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where  i  and  j  are  indices  for  the  matrix  roles  and  columns,  respec- 
tively,  (1,1  being  at  the  top  left  part  of  the  matrix)  and  where  A*  is 
just  the  matrix  A  multiplied  by  itself  x  times.  The  coefficients  of  A* 
have  well-defined  structures,  easily  obtained  by  simply  carrying  out 
the  indicated  multiplications.  For  those  cases  where  the  number  of 
launch  attempts,  x,  is  larger  than  r,  where  r  is  the  rank  of  the  matrix 
A,  the  first  r  -  1  coefficients  of  the  first  column  in  A*  are  simply  the 
first  r  -  1  terms  in  the  expansion  of 

(Qs  P3)  ' 


respectively.  The  rth  term  is  simply  one  minus  the  sum  of  the  first 
r  -  1  terms.^  The  coefficients  of  the  second  column  are  the  same  as 
the  first,  slid  down  one  row,  with  the  last  coefficient  equal  to  one 
minus  the  sum  of  the  first  r  -  2  terms,  and  so  forth  for  the  remaining 
columns.  Of  course,  the  last  column  is  all  zero*;  \cept  for  the 
coefficient  at  the  bottom. 

Let  j(x)  be  the  coefficients  of  the  matrix  A* .  Then  assuming  that  x 
>  r,  the  coefficients  are  as  follows: 

ai.|(x)  =  ISi<r,  lsj<i  (A.81) 


where  .  ^ .  is  the  binomial  coefficient: 
\*  K 


aj.j(x)  =  0 

for  all  i  <  j  <  r: 

(A.82) 

r-l 

ar.jW  =  l-£aj.j(x) 

'=i 

for  ail  j  <  r: 

(A.83) 

and 

r(x)  =  I 

(A.84) 

•’The  sum  of  all  the  coefficients  in  each  column  of  the  matrix  must  be  equal  to  one. 
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The  coefficients  for  x  <  r  are  equally  straightforward,  but  we  will  not 
offer  them  here.  Figure  A.1  shows  the  domains  of  the  matrix  associ¬ 
ated  with  each  of  the  above  four  equations. 

The  utility  of  the  above  solution  is  questionable,  and  calculations 
have  simply  used  Eqs.  (A.74)  through  (A.78). 

SECTION  THREE:  THE  MONTE  CARLO  APPROACH  TO 
SOLVING  THIS  PROBLEM 

The  third  approach  that  can  be  taken  is  to  simply  write  a  Monte 
Carlo  computer  program  where  each  TEL  is  subjected  to  prelaunch 
and  posdaunch  counterforce  attacks.  The  flow  chart  of  the  program 
is  shown  in  Figure  A.2. 

On  each  trial,  the  available  TELs,  loaded  with  TBMs,  sortie  from  their 
hide  sites  to  their  launch  sites.  Random  numbers  are  drawn  for  each 


.  RAN01UA469-A  f 

Matrix  A* 


Figure  A.1 — Domains  of  Matrix  A  Associated  with 
Eqs.  (A80)  Through  (A.83) 
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Figure  A.2 — Flow  Chart  for  Monte  Carlo  Calculation 


TEL  to  determine  whether  the  prelaunch  counterforce  attack  against 
that  TEL  succeeds.  Assuming  that  the  TEL  survives,  it  launches  its 
TBM  unmolested.  Subsequent  to  the  launch,  the  TEL  is  subjected  to 
a  counterbattery  attack.  A  new  random  number  is  drawn  for  each 
TEL,  determining  whether  it  survives  that  attack.  Those  TELs  that 
survive  return  to  their  hide/resupply  site,  reload,  and  are  available  to 
participate  in  the  next  sortie.  The  numbers  of  TBMs  launched  are 
finished  when  one  of  three  conditions  is  met:  (1)  all  available  TBMs 
are  either  launched  or  killed  during  the  prelaunch  attack  phase,  (2) 
all  TELs  are  dead,  or  (3)  the  maximum  number  of  trials  (specified  by 
the  user)  has  been  reached  (not  shown  on  Figure  A.2). 

Multiple  trials  are  run,  usually  somewhere  between  100  and  1000 
separate  runs,  to  obtain  average  outcomes  not  significantly  different 
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from  the  expected  values.  In  addition  to  the  statistical  averages,  the 
standard  deviation  of  the  outcomes  horn  the  mean  is  also  calculated. 
Table  A.1  shows  an  example  of  the  data  obtained  by  running  the 
model.  It  shows  that  after  500  trials,  the  model  estimated  that  on  the 
average  49.57  TBMs  would  be  launched,  given  the  input  assump- 

TableA.1 

Sample  Outcome  from  Monte  Carlo  Run  for  Counterbattery  Attack 


Input  Data 

Total  TELs  in  cluster  =  100 
Total  TBMs  in  cluster  =  10 
Max.  salvos  =  50 


Pfcvs.  TEL  =  0.2 
Number  of  trials  =  500 

Output  Data 

Average  number  of  TBMs  fired  =  49.57 
Standard  deviation  =  13.241 
Standard  deviation  mean  =  0.5922 


Histogram  of  Number  of  TBMs  Fired 
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lions.  This  outcome  is  close  to  the  value  obtained  by  other  analytic 
means.  The  estimated  deviation  of  the  results  about  this  means  is 
estimated  to  be  13.241  launches,  a  substantial  number.  Note  that  the 
model  also  estimates  the  standard  deviation  of  the  estimate  of  the 
mean — how  far  off  our  estimate  is  likely  to  be. 

The  histogram  of  the  number  of  TBMs  launched  is  also  shown.  The 
first  column  represents  the  number  at  the  start  of  the  row,  and  the 
numbers  in  the  top  row  represent  the  number  to  be  added  to  the  first 
column.  Thus,  to  find  out  how  many  times  exactly  55  TBMs  were 
launched,  read  down  the  first  column  to  50,  then  across  the  first  row 
to  5.  You  should  find  that  for  these  500  trials  exactly  13  resulted  in 
producing  55  TBM  launches.  It  is  interesting  to  note  that  one  trial 
produced  100  launches,  a  very  unlikely  outcome,  and  the  minimum 
number  of  launches  was  18,  also  a  very  unlikely  outcome  given  that 
the  first  10  launches  occur  before  any  counterbattery  attacks  occur.*’ 


®For  these  input  assumptions,  there  is  less  than  one  chance  out  of  four  that  in  500 
trials  as  few  as  1 8  TBMs  would  be  all  that  were  launched. 


_ Appendix  B 

COMPUTER  PROGRAMS  FOR  CALCULATING 
OUTCOMES  FROM  PRELAUNCH  AND  POSTLAUNCH 
COUNTERFORCE  ATTACKS  AGAINST  TBMs  AND  TELs 


1.  Simple  counterforce/counterbattery  model  to  calculate  TBMs 
launched  in  the  single  TEL  case. 

CLS 

OPEN  "cUp;"  FOR  OUTPUT  AS  #1 

'input  pre-launch  Pk,  vary  over  counterbattery  Pk 

INPUT  "Prob  of  pre-launch  counterforce  kill  (PI)  ";pl 
DIM  E(20) 


'start  calculation  at  Pk  less  than  one  to  avoid  zero  divisor 
p2=.9 


WHILE  p2>0 


ql=l-pl 

q2=:l-p2 

q3=ql‘'q2 


'Calculate  for  1-10  TBMs  per  TEL 
FOR  N=1  TO  10 

E(N)=ql*(l-q3''N)/(l-q3) 

NEXTN 

E(ll)=l/(l-q3) 


'Write  to  clipboard  for  later  inclusion  in  EXCEL  or  CRICKET  graphs 

WRITE  #1,  pl,p2,E(l),E(2),E(3),E(4),E(5),E(6),E(7),E(8),E(9),E(10),E(ll) 
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'Adjust  step  size  to  reflect  sensitivities  at  small  Pk 

k=.l 

IF  p2<.3  THEN  k=.02 
IF  p2<.08  THEN  =.01 
p2=p2-k 

WEND 

CLOSE  #1 

2.  BASIC  program  to  calculate  expected  number  of  TBMs 
launched,  given  a  postlaunch  counterbattery  attack. 

This  program  does  the  following  things: 

1.  It  is  based  on  mathematical  equations  for  counterbatterv  attacks 
only. 

2.  Intermediate  steps  calculated  the  probabilities  that  specific 
outcomes  are  achieved.  This  information  is  not  available  in  other 
calculations. 

3.  The  program  is  limited  to  six  TELs.  Generalizations  beyond  six 
are  straightforward,  but  the  compute  time  would  grow 
significantly. 

4.  Prints  results  to  clipboard,  to  permit  use  of  EXCEL  or  CRICKET 
graphics. 

5  Automatically  calculates  outcomes  for  various  Pks. 

'  Program  to  Calculate  Expected  Value  Outcome  of  Counterbattery 
Attacks 

CLS 

OPEN  "CLIP: "  FOR  OUTPUT  AS  #1 
DIM  zl(200) 

DIM  z2(200) 

DIM  z3(200) 

DIM  z4(200) 

DIM  7.5(200) 

DIM  z6(200) 

DIM  z7(200) 

DIM  z8(200) 
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'Input  variables 

INPUT  "Number  of  TBMs;",  N 
x=0 

10  Pk=.9 

'Start  of  calculation 
WHILE  Pk>0 
x=x+l 


'Calculation  for  1  TEL 

El=(l-(l-Pk)^N)/Pk 


'Calculation  for  2  TELs 

B2=0 

FOR  K=2  TO  N-1 
B2=B2+(N-K)*(K-  l)*(l-Pk)''(K-2) 
NEXTK 
E2=N-B2*Pk'^2 

'Calculation  for  3  TELs 

B3=0 

FOR  K=3  TO  N-1 

B3=B3+(N-K)’'(K-l)»(K-2)»(l-Pk)^(K-3) 

NEXTK 

E3=N-(B3*PkA3)/2 

'Calculation  for  4  TELs 

B4=0 

FOR  K=4  TO  N-1 
C4=0 

FOR  i=l  TO  K-3 
C4=C4+(K-l-i)’'(K-2-i) 

NEXTi 

B4=B4+(N-K)*C4*(l-Pk)^(K-4) 

NEXTK 

E4=N-(B4*Pk''4)/2 


'Calculation  for  5  TELs 

B5=0 

FOR  K=5  TO  N-1 
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C5=0 

FOR  i=l  TO  K-4 
D5=0 

FORj=lTOK-3-i 

D5=D5-E(K-l-i-j)*(K-2-i-j) 

NEXTj 

C5=C5-hD5 

NEXTi 

B5=B5-i-(N-K)»C5*(l-Pk)''(K-5) 

NEXTK 

E5=N-(B5*Pk^5)/2 


'Calculation  for  6  TELs 

B6=0 

FOR  K=6  TON-1 
C6=0 

FOR  il=l  TO  K-5 
D6=0 

FOR  i2=lTOK-4-il 
E6=0 

FOR  i3=l  TO  K-3-il-i2 
E6=E6+(K-l-il-i2-i3)»(K-2-il-i2-i3) 
NEXT  i3 
D6=D6-t-E6 
NEXTi2 
C6=C6-^D6 
NEXTil 

B6=B6-i-(N-K)»C6»(l-Pk)''(K-6) 

NEXTK 

E6=N-(B6»Pk^6)/2 


'Output  routines,  print  both  to  screen  and  to  clipboard 

PRINT  USING  "###";  N, 

PRINT  USING  "##.##■•;  Pk, 

PRINT  USING  "###.###";  E1,E2,E3,E4,E5,E6, 

PRINT  USING  "##.##";  EUPk,  E2''Pk,  E3»Pk,  E4*Pk,  ES^Pk,  E6»Pk 

zl(x)=N 

z2(x)=Pk 

z3(x)=El 

z4(x)=E2 


I 
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25{x)=E3 

z6(x)=E4 

z7(x)=E5 

z8(x)=E6 

WRITE  #1,  zl{x  ^4(x),  z5(x),  z6(x),  z7(\),  z8(x) 

'change  value  of  Pk  and  repeat  calculation 
IF  Pk>.2  THEN  del=.l  ELSE  del=.02 
IF  Pk<.06  THEN  del=.01 
Pk=Pk-del 

WEND 

'Request  new  inputs  to  run  another  case 
INPUT  "change  number  of  TBMs";  s$ 

IF  s$="y"  GOTO  100 
GOTO  110 

100  INPUT  "Number  of  TBMs:",  N 
GOTO  10 

'End  program 
110  CLOSE  #1 

3.  Markov  model  program  to  obtain  expected-value  outcomes  for 
prelaunch  counterforce  and  postlaunch  counterbattery  attacks 


CLS 

OPEN  "clip:"  FOR  OUTPUT  AS  #1 

'Input  values  for  calculation 

INPUT  "Number  of  TELs=";  NuTEL 

INPUT  "prelaunch  prob  of  kill=";pl 

INPUT  "post-launch  prob  of  KILL=";p2 

INPUT  "Maximum  number  of  TBMs=";N 

q=(l-pir(l-p2) 

p=l-q 

'Dimension  arrays 
DIM  T(20,100) 


I 
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DIM  L(IOO) 
DIM  Tot(lOO) 


'Set  initial  conditions,  T  is  prob  state  of  TEL  set 
FOR  i=OTONuTEL-l 
T(i,0)=0 
NEXTi 
T(NuTEL,0)=l 
x=0 


WHILE  x<N:  'Calculate  outcome  for  all  x<=N 

X=X-El 

L(x)=(l-T(0,x-l))*(l-pl):  'L  is  probability  of  a  TBM  launch  for  each  x 
Tot(x)=Tot(x-l)-EL(x):  Tot  is  cumulative  number  of  TBMs  launched 
at  X 

'Matrix  calculation  for  TEL  prob  state 
T(NuTEL,x)=qT(NuTEL,x-l) 

FOR  i=l  TO  NuTEL-1 
T(i,x)=q*T(i,x-l  )+p'*T(i-fl,x-l ) 

NEXTi 

T(0,x)=T(0,x-l)-HpT(l,x-l) 

'Show  outcomes  for  each  x 
PRINT  X,  L(x),  Tot(x),  T(0,x), 

WRITE  #1,  x,L(x),Tot(x) 

WEND 

CLOSE  #1 

4.  Monte  Carlo  program  for  assessing  outcomes  of  combined 
prelaunch  counterforce  and  postlaunch  counterbattery  attacks 
on  reducing  the  total  number  of  TBMs  launched. 

'Counterforce  and  Counterbattery,  MC.4 


I 
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This  progam  does  the  following  things: 

'  1.  In  monte  carlo  mode  calculates  number  of  TBMs  launched, 

’  assuming  pre-  and  post-launch  counterforce  attacks  by  the 
defended. 

2.  Runs  in  automatic  mode  500  trials  for  each  case. 

3.  Performs  cases  for  input  values  for  TBM  and  TEL  inventories. 

4.  Performs  cases  for  input  values  for  pre-  and  post-launch 
counterforce  effectiveness 

(Pkl  and  Pk2) 

5.  Calculates  expected  outcomes  and  their  standard  deviations. 

6.  Keeps  statistics  on  outcomes  per  trial,  prints  distributions  of 
outcomes. 

'Start  of  program 

CLS 

INPUT  "Specify  name  of  output  file z$ 

OPEN  z$  FOR  OUTPUT  AS  #1 
DIM  T(IOOO) 

DIM  Attack(lOOO) 

DIM  Outcome(lOOO) 

DIM  c(lOOO) 

DIM  d(10(K)) 

DIM  liveTEL(lOOO) 

RANEXDMIZE  TIMER 

'Inputs  for  case 

INPUT  "Total  TBM  inventory  per  cluster";  totTBM 
INPUT  "Total  number  of  TELs  per  cluster";  totTEL 
INPUT  "Pre-launch  Counterforce  Effectiveness";  Pkl 
INPUT  "Post-launch  Counterbattery  Effectiveness";  Pk2 
INPUT  "Number  of  random  trials  desired,  <Retum>  sets  no.  at  500"; 
Tr 

IF  Tr=0  THEN  Tr=500 

'Prints  (to  readable  file)  input  variables  at  the  start  of  each  case 
PRINT  #1,"  " 

PRINT  #1,  "Total  TELs,  TBMs  in  a  cluster";  totTEL;  totTBM, 

PRINT  #1,  "Max  Salvos";  W;  "Pkl";  Pk,"Pk2",Pk2,  "Num  of  trials";Tr 


I 


I 
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'Initializes  variables  at  the  start  of  each  case 
FOR  y=0TO999 
c(y)=0 
d(y)=0 
NEXTy 
x=0 

umax=0 

SuinAttack=0 

SumVar=0 

'Performs  calculations  for  TBMs  launched  for  Tr  trials 
FOR  x=l  TO  Tr 

'Set  up  initial  conditions  for  calculations  of  each  trial 

liveTEL(x)=totTEL 

NuTBM=totTBM 

Attack(x)=0 

'Checks  to  see  whether  calculation  should  stop 
WHILE  NuTBM*UveTEL(x)>0 

'Pre-launch  attack  calculation 

NuTBM=NuTBM-l 

a=RND 

IF  a<Pkl  THEN 
liveTEL(x)=IiveTEL(x)-l 
GOTO  70 
END  IF 

'Accumulates  total  TBM  launches 
Attack(x)= Attac  k(x)-El 

'Post-launch  attack  calculation 
a=RND 

IF  a<Pk2  THEN  liveTEL(x)=liveTEL(x)-l 
70  WEND 


I 
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'Summary  data  from  trial 

SumAttack=SumAttack+Attack(x) 

SumVar=SumVar+Attack(x)''2 

'Prepares  histogram  for  total  TBMs  launched  per  trial 

y=Attack(x) 

c(y)=c(y)+l 

'Prepares  histogram  for  surviving  TELs 

z=liveTEL(x) 

d(z)=d(z)+l 

120  NEXT  X 

'Calculations  for  trials  can  take  awhile,  alerts  user  case  done 
BEEP 

'Prepares  summary  and  distributional  output  data  for  each  case 

Ave=SumAttack/Tr 

Var=(SumVar/Tr)-Ave''^2 

Sigma=SQR(Var*Tr/(Tr-l)) 

Sig.mean=Sigma/SQR(Tr) 

'Prints  data  to  readable  file 

PRINT  #1,  "Trials=";Tr;’Ave  No.  TBMs  fired=";Ave,  "Std 
Dev=";Sigma,"Std  Dev  Mean=";Sig.mean 
PRINT  #1,"  " 

PRINT  #1,  "Histogram  of  number  of  TBMs  fired" 

PRINT  tfl, " 

PKir  JT  #1,  USING  "#####",-0;1A3;4;5;6;7;8;9 
''OR)=  TO  totTBM  STEP  10 

fi<(N  1  #1,  USING  "#####",-j;c(j);c(j+l);c(j+2);c0+3);c(j+4);c(j+5); 
c(j+6);c(j+7);c(j+8):c(j+9) 

NEXTj 
PRINT  #1,"  " 

PRINT  #1,  "Histogram  of  number  of  TELs  sur\'iving" 

PRINT  #1," 
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PRINT  #1,  USING  "#####"A1,2,3A5,6,7,8,9 
FOR  j=0  TO  totTEL  STEP  10 

PRINT  #l,USING"#####’*;j;d{j);d(j+l):d(j+2);d(j+3);d(j+4);d(j-t-5); 
d{j+6):d(j+7):d(j+8):d(j+9) 

NEXT) 

'Closes  the  program 

CLOSE  #1 

END 


